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ABSTRACT 

We confirm our earlier tentative detection of M31* in X-rays and measure its 
light-curve and spectrum. Observations in 2004-2005 find M31* rather quiescent 
in the X-ray and radio. However, X-ray observations in 2006-2007 and radio 
observations in 2002 show M31* to be highly variable at times. A separate 
variable X-ray source is found near PI, the brighter of the two optical nuclei. 
The apparent angular Bondi radius of M31* is the largest of any black hole, and 
large enough to be well resolved with Chandra. The diffuse emission within this 
Bondi radius is found to have an X-ray temperature ~ 0.3 keV and density 0.1 
cm -3 , indistinguishable from the hot gas in the surrounding regions of the bulge 
given the statistics allowed by the current observations. The X-ray source at the 
location of M31* is consistent with a point source and a power law spectrum 
with energy slope 0.9 ± 0.2. Our identification of this X-ray source with M31* is 
based solely on positional coincidence. 
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Introduction 



SMBHs in galactic nuclei spend the vast majority of their life accreting at very low 
rates, but our understanding of how this accretion occurs is relatively poor compared to our 
understanding of what happens at high rates. At high rates a slim accretion disk often forms 
and ~10% of the accretion energy is radiated. At low rates the accretion becomes radiatively 
inefficient, but we are uncertain what fraction of the accretion energy is radiated, and what 
fraction of the gas accreted at large radii actually reaches the black hole event horizon. 

While it is clear that some sort of a radiative ly inefficient accret io n flow (RIAF) oc- 
curs, its form as a magnetically dominated in flow IShvartsman (119711 ): iMelial (1l992aJ ) . an 
ADAF (Advection Dominated Accretion Flow; Narayan fc Yi ( 19941 )). CDAF (Convective 
Dominated Accretion Flo w in which ady e ction occurs but is moderated by convection; 
Igumenshchev et all (boooh : iNaravan et al.l (hood ): buataert fc Gruzinovi (bood )). an ADAF 
& wind (in which a large fraction of the accreted material is blow out in a wind before 
reaching the SMBH; Narayan fc Yi J 19951 )); sometimes c alled an ADIOS (Advecti v e Dom - 
inated Inflow Outflow Solution; iBlandford fc Begelmanl (119991 ); lHawley fc Balbud (120021 )). 
or something else is unclear. One of the most promising ways to determine which of these 
alternatives is correct is to image the accretion flow on a sub-Bondi scale and therefore de- 
termine the structure of the flow. Spectra taken on sub-Bondi scale would also be highly 
discriminating. 

Perhaps the best candidate for such a study is the SMBH in M31, called hereafter M31* . 
The mass of M31* has previously been estimated as 3 x 10 7 M^ dKormendv &: Bender 1999), 



but more recent HST spectroscopy indicates the mass is lAtol x 10 8 M^ dBender et al.l 



2005). 



This higher mass means the Bondi radius is not the 0.9" previously estimated (IGarcia et al. 
20051 ) but ~ 5" (see below), making it the largest of any known SMBH (see Figure 1). 



Within the Bondi radius Chandra det ects several point sources and also diffuse gas which is 
presen t throughout the bu lge of M31 (Shirev et al.ll2001at iTakahashi et al.ll2004l ; iLi fc Wang 



2007al ). Our earlier work (IGarcia et al.ll2000l ) suggested that a super-soft point source near 
the nucle us was the X-ray counterp art of M31*, but subsequent alignment with M31 globular 



clust ers (IBarmby &: Huchral l200ll ) revealed that this was not the case ( IGarcia et al.l 12001 



20051 ). Chandra HRC imaging revealed that the X-ray source closest to the nucleus is made 
up of two partially resolved components, one of which is within ~ 0.1" of the position of 



P3 and which we t entatively identified as the ~ 10 36 erg s 1 X-ray counterpart of M31 



(IGarcia et al.ll2005l ). In this work we confirm the detection and measure the variability of 
this SMBH. 



Hubble imaging revealed an unusual double nucleus within M31 (ILauer et al.l Il993l ) 
which has been successfully modeled as an eccentric torus of stars viewed nearly edge on 
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( ITremaindll995l ) . The optically brighter end of the torus is known as PI, the optically fainter 
(but UV brighter) end as P2. More recent Hubble imaging and spectroscopy revealed that 
embedded within P2 is a dense cluster of several hundred A-stars which are very UV bright. 
This cluste r has been dubbed P3 and its brightness peaks at essentially the same position as 
that of P2 (IBender et al.ll2005l ). HST spectroscopy reveals a dark mass (presum ed black hole) 
located at the center of this cluster with an estimated mass of 1.4±g;| x 10 s M^ flBender et all 



20051 ). The presence of so many young A stars in the nucleus of M31, which otherwise is 
made u p of an old stellar p opulation, is a myste r y. Th e same phenomenon is also seen near 
Sgr A* (jGhez et~aL 2003a ). Demarque fc Virani ( 2007 ) suggest that these apparently young 
stars may actually be old and relics of stellar collisions in the very dense nuclear regions, but 
recent HST spectra in dicate that at least in the case of the Ga lactic center these are genuine 
young, massive stars (jGhez et al.l l2003bl ; iMartins et al.l 120081 ). A suggested source for the 
gas which later collap ses to form the P3 star cluster is stellar mass loss from the older stars 
in the orbiting torus (jChang et al.l 120071 ). 



Spitzer imaging reveals that M31 looks less like a spiral galaxy than previously thought. 
Instead, the warm gas and dust is in the form of a set o f rings, apparently due to transits 



of M 32 through the center of the plane of the M31 disk ( Barmby et al.l 12006; 



Gordon et al. 



200 



3), 



20061 ). The most recent of these transits took place 210 million years ago (IBlock et al 
and would have triggered a burst of star formation. 

The associatio n of the ra d io sour ce (~ 30 /iJy at 3.6 cm) with a nuclear black hole by 



Crane et al.l (119931 ) and iMelial (Il992bl ) and subsequently by all others so far is supported by 
the detection of an unre solved (< 1 pc or .35") radio source at the position of the stellar 
nucleus (within ~ 1 /2". iGarcia et al.l (120051)). and radi o variability which is comparable to 



the radio variability of Sgr A* (jSjouwerman et al.ll2005l ). 



While M31* is lOOx further away than the SMB H in the center of our Galaxy, it suffers 

1 JGarcia et all l200oh . whereas A v ~ 3 for Sgr A * Th e 



much less reddening: A 



v 



Shirev et al. 


( 


2001b 


); 


Li & Wane f 


2007b): 



.3 keV (iDosai et all (l200lf ): 



(M31* and Sgr A*) will deepen our understanding of the physics of black hole accretion 
and emission mechanisms at low rates. In this regard, M31* and Sgr A* offer interesting 
differences in that M31* is 3x less luminous in the radio but 1 to 3 orders of magnitude more 
luminous in the X-ray 
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Fig. 1. — The Bondi radii of nearby SMBH vs. their apparent X-ray luminosity (or upper 
limits) in units of the expected Bondi luminosity. Objects in the lower right of this diagram 
provide simultaneously secure and severe constraints on accretion and emission models for 
these SMBH. In this regard, M31* is the outstanding object. 
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Observations 



M 31 has been observed e x tensively with both Cha ndra (i.e.. lGarcia et al.l (120001) ;lKong et al 



(2003 


); 


Williams et al. 


(2001b) 


; Pietsch et al. 



(2006): iDi Stefano et al.l (120041 )) and XMM- Newton 



i.e.. 



Shirey et al. 



(120071 )). Many of the Chandra observations have been in the form of 



monthly short (~5 ks) snap-shots in order to monitor the variability of the point sour ce pop 



ulatio n, but there also have been longer observations to search f o r rapid variability (Kaaret 



20021 ) and to study source populations (IDi Stefano et al.l l2004t iPietsch et al.l 120071 ). The 
total Chandra exposure on M31 is now ~ 1.3Ms and it continues to climb. If one limits 
the selection of data to only those taken with ACIS-I within 1 arcmin of the nucleus, the 
total exposure is currently ~ 250 ks. While Chandra is able to view M31 10 months out 
of every year, visibility from XMM- Newton is restricted to a few week long window every 6 
months, so the XMM- Newton obser vations have taken the form of less frequent, but longer, 
exposures (jTrudolyubov et al.ll2005l ). 



The Chandra observations discovered the first resolved SNR in an external galaxy 
(IKong et al.ll2002al ). and summed Chandra observations have provided a fiducial point for 
X-ray color-color diagrams which help to separate unresolved X-ray source populations in 
more distant galaxies into SNR, high and low mass X-ray binaries, and background AGN 
(IKong et al.l l2002bl ; iPrestwich et al.l 120031) . These da t a-sets have also revealed extensive 
structure in the diffuse emission (ILi fc Wangj l2007bl ; iBogdan fc Gilfanovl 120081 ). Multi- 
wavelength studies of the M31* environm ent suggests the presence of an outflow of hot 
gas from the nuclear region ILi et al.l (120091 ) . 



Below we discuss the data-sets analyzed in this paper; first those taken with the Chandra 
HRC and then those taken with the Chandra ACIS. 



2.1. Chandra HRC Observations 



In 2004/2005 we undertook a series of 4 moderate (50ks) HRC-I exposures separated 
by 1 month in order to confirm our earlier possible detection of M31* and to search for any 
variability. Two of these exposures were subdivided into approximately equal sections and 
separated by the normal 10 hour pause in observations while Chandra transits the Earth's 
radiation belts. We also include in our analysis archival HRC observation s made during 
2006/2007 which were design ed to monitor optical nova (IHenze et al.l 120091 ). and a single 
50 ks observation from 2001 (IKaaretl 120021 ) . Table 1 lists the data-set we have analyzed, 
giving observation dates and exposure times. 



In order to identify M31* within the Chandra images and then extract a light-curve, it 
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Table 1: This table shows the Observation ID's, Dates and Exposure Times for the 16 
Chandra HRC images discussed in this paper. 



was first necessary to accurately determine the location of M31* within the HRC images. 
The absolute astrometric accuracy of HST and Chandra images is < 1", which is insufficient 
to uniquely identify sources in the crowded core of M31. Within n arrow fields ( 1 1') it i s 
possible to register HST and Chandra images to ~ 0.1" accuracy (jEdmonds et al.ll2003l ). 
Over larger fields the PSF of the Chandra mirror increases substantially, and this along 
with the calibration of the HRC may limit the astrometric accuracy even with registration 
( IBeckerman et al.ll2004l ). 

There are very few sources in common between our Chandra images and the HST/ACS 
(F435W) images which show the M31 double nucleus, so we used a 2 step process to reg- 
i ster our images. As part of our search for counterparts to black hole X-ray transients 
( jWilliams et al.ll2005l ) we obtained several HST/ACS i mages of the M31 n ucleus, and one 
of these was first registered to the Local Group Survey (IMassey et al.ll2006l ) image of M31. 
The second step is to register the Chandra images to the same LGS image. We allow for an 
x and y translation, a single change in scale factor, and roll, but typically the latter two of 
these are negligible while the translation is significant and < 1". 
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The HST/ACS to LGS registration typically uses 20 stars and has a RMS of 0.03", and 
the registration of individual Chandra/HRC images to the LGS uses 6 or 7 globular clusters 
and has a rms of 0.1" — 0.2". One of these 7 globular clusters sometimes falls below the 
Chandra detection threshold, limiting us to 6 occasionally. We note that there are 9 globular 
clusters which are in both the Chandra and LGS images, but 2 of these are between 6' and 
10' from the nucleus where the PSF due to the Chandra mirror has widened substantially. 
Including these distant clusters increases the error in the registration to ~ 0.3", so we have 
excluded them. The remaining clusters are between 1' and 5' from the nucleus. Clearly the 
Chandra to LGS registration error dominates the final registration, and this error itself is 
dominated by the counting statistics within the PSF (and associated registration errors) in 
the individual globular clusters within the Chandra images. 

In order to ensure the most accurate, and perhaps more importantly, stable, registration, 
we summed the HRC data into a single image and registered that image to the LGS images. 
The summing process (CIAO reproject-aspect) uses the positions of bright X-ray sources to 
register the individual HRC images to a common frame to better than 0.05". The accuracy 
of the registration of the summed Chandra image to the LGS image as determined with 
IRAF ccmap is 0.1" rms. We also summed the 2001/2004/2005 and 2006/2007 HRC data 
separately and registered these two images to the LGS, again finding a rms accuracy of 0.1". 
Comparing these two (now registered) HRC images to each other, we find an offset between 
bright sources of 0.2", larger than the formal error in the rms. Given that we are solving for 
4 parameters (x and y shift, scale factor, and roll) with 7 points, the resulting formal rms 
might not be expected to be Gaussian - we therefore take 0.2" as our de-facto registration 
error between the HST and Chandra images. 

These two images are shown in Figure 2. The 2006/2007 image shows two separate X- 
ray sources at the approximate locations of PI and P3. A simple translation of 0.25" brings 
the Southwestern source to the position of P3, and the Northeastern source 0.1" to the South 
of PI. Given the concurrence between the separation of the X-ray sources, the angle between 
them, and PI and P3 as seen in the HST images, we identify the X-ray sources with M31* 
(=P3) and a source within PI. This concurrence can also be seen in Figure 3, which shows 
the HST/ACS image with the HRC contours overlaid. However, we cannot entirely exclude 
the possibility that one or more transient sources unrelated to PI or P3 are responsible for 
the X-ray emission. 

Since counting statistics of the globular clusters in the individual ~ 20 ks X-ray images 
limit the accuracy of the registration, we rely on the registration of the merged 430 ks image 
and measure offsets to PI and P3 from the nearby super-soft source (CXO J004244. 2+411608, 



Garcia et al.l (120001 ); hereafter SSS). The SSS has between 100 and 300 counts in each of the 
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Fig. 2.— The summed 2001/2004/2005 (left) and 2006/2007 (right) HRC images of the M31 
nucleus, registered to a common astrometric frame. The registration is accurate to ~ 0.2". 
The contours are from an HST/ACS image. The innermost contour levels are closed at the 
positions of PI (upper left, or North-East) and P3 (lower right) and are separated by 0.5". In 
the 2004/2005 image the source associated with P3 (=M31*) is faint, while in the 2006/2007 
image this source is stronger and clearly separate fr om the source near PI. The super-soft 
source CXO J004244.2+411608 JGarcia et all boooh is the bright source directly South of 
P1/P3. 



images, so can be centroided to ~ PSF/ ^/{counts) ~ 0.05" in each of the individual images. 

In order to extract the HRC light-curves we centered extraction circles at the locations 
of PI and P3 as determined by the offsets from the SSS (after the translation of 0.25" 
discussed above) on each of the images and counted the photons therein. We used radii of 
0.4" centered at PI and P3, as shown in Figure 4. PI and P3 are sufficiently close that these 
circles overlap, so we excluded the counts in the overlapping region. 

Figure 5 shows the resulting light-curves. There is clearly some correlation between the 
PI and P3 light-curves, as is expected due to the proximity of the two sources. The 0.4" 
extraction circle contains 50% of the flux from a point source and the two sources are only 
0.5" apart. 

The most dramatic variability we find is in December 2007 (day 519 to 529), where 
M31* varies by a factor of 3 in 10 days. The maximum variability is between day 117 and 
519, where we see a factor of more than 10 variation. In order to search for variability 
on short timescales, we divided the day 519 observation (where M31* was brightest) into 5 
approximately equal 4000 s intervals. Figure 6 shows the resulting light-curve. While there 
are apparent ~ 2a deviations, a \ 2 test to the hypothesis of a constant source has a 6% 
chance of being satisfied, so significant variability cannot be convincingly argued. 
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Fig. 3. — A false color representation of the HST/ACS F435W image of the M31 nucleus, 
with the contours from the HRC images overlaid. The left image shows the contours from 
the summed 430ks HRC data. The right image shows the contours from the 2006/2007 HRC 
data alone. A small East- West translation will bring the X-ray peaks on top of P3 and 
within 0.1" of PI. 

In order to test the sensitivity of the light-curve to possible errors in the registration, 
we shifted the extraction circles 0.05" both RA and DEC, as this is the typical error in 
measuring the centroid of the SSS which is used as the reference point. While the overall 
character of the resulting light-curve remains the same, some of the points do see > 2a 
changes (see Figures 5 and 6). 

Assuming the spectral shape found with ACIS (see below) of a u = 0.9 and N H = 6 x 10 20 
cm~ 2 , a measured rate of 1.0 c/ks corresponds to an emitted luminosity over 0.3 — 7.0 keV 
of 2.25 x 10 36 erg s" 1 in the full HRC beam, or 3.0 x 10 36 erg s" 1 over 0.1 - 7.0 keV. M31* 
therefore ranges in luminosity from 0.6 — 20 x 10 36 erg s _1 (0.1 — 7.0 keV). We note that if 
the source we associate with M31* is in fact unrelated, then these number are upper limits 
and the X-ray luminosity of M31* must at times be below 6 x 10 35 erg s _1 . 



Fig. 4. — Two of the 15 HRC images used in our light-curves. The left image from 2004 
Dec 6.75 shows the source at PI somewhat brighter than that at P3, while the right image 
from 2008 Dec 7 shows P3 at its brightest. The 0.4" circles used to measure the counts are 
re-centered for each observation by using an identical set of vect ors positioned to o riginate 
at the centroid of the super-soft source CXO J004244.2+411608 jOarcia et al.lboool ) seen at 
the bottom of the images. Light curves were generated from the counts in these 0.4" circles, 
excluding counts from the overlapping region. 
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Fig. 5.— Light-curves for P3 = M31* and the source near PI as extracted from our 
2004/2005 HRC data (left) and archival 2006/2007 (right) HRC data. A counting rate of 1.0 
counts per ks corresponds to a luminosity of 3 x 10 36 erg s _1 (0.1 — 7.0 keV) at the 780 kpc 
distance of M31. Error-bars are la counting statistics only. The dashed line indicates the 
sensitivity of the light-curve to la shifts in the registration and therefore count extraction 
procedure. We do not include the 2001 observation in the light-curve due to its 3 year offset. 
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Fig. 6.— Light-curve for P3 = M31* during obsid 8529 (day 529 in Figure 5). The data 
were divided into 4000 sec intervals in order to measure the rate. While there are ~ 2a 
deviations, a \ 2 test to a constant source finds y^/v = 8.87/4, which has a 6% chance of 
originating from a constant source. Short time scale variability is therefore not detected. As 
before, the dashed line and points represent the light-curve extracted from slightly shifted 
positions. 
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2.2. Chandra ACIS 

While the HRC data gives excellent spatial resolution and allows us to determine an 
accurate light-curve for M31*, They do not provide sufficient energy resolution to allow us to 
search for spectral signatures of the Bondi flow/ADAF. We therefore investigate the archival 
AC1S-I and AC1S-S data. 

We used 29 ACIS-I exposures aimed within 1' of M31* which total to 148ks of exposure 
time. This summed exposure is shown in Figure 7 (left) as a smoothed 3-color X-ray image 
where the soft (0.3-1.0 keV), medium (1.0-2.5 keV), and hard (2.5-7.0 keV) bands are color 
coded as red, green, and blue respectively. The 5" Bondi radius is shown here in green. The 
emission due to hot diffuse gas is seen throughout this image and generates ~ 1 event per 
pixel in this image. Note that we do not include the ACIS-S observation discussed below, 
obsid 1575, in this sum due to the substantially different energy response of ACIS-S vs. 
ACIS-I. 

We searched for an enhancement in the flux at (or inside) of the Bondi radius by 
generating radial profiles of the flux, centered on M31* , in several wedges that are free of 
point sources. Defining the ray to the North as zero degrees and moving counterclockwise, 
we generated these profiles over the azimuth ranges of 40 to 70 degrees, 100 to 150 degrees, 
and 245 to 280 degrees. In the first two of these azimuth ranges we found a slight excess 
of counts near the Bondi radius (above the local background), amounting to 10 counts from 
3" to 5" in the first region and 20 counts from 4" to 6" in the second region. We then 
used MARX to estimate the number of counts that would be scattered into these regions 
from the surrounding point sources, and found that this could account for ~ 1/3 of the 
observed excess. No te that there is clear spatial structure in the diffuse emission at larger 



radii (10" — 20", see iLi et al.l (120091 )) which, if present at or within the Bondi radius, could 
easily account for the apparent excess of 10 to 20 counts seen in two of our regions. Assuming 
that the results in these three azimuth ranges are typical of what would be found for the 
diffuse emission at all azimuth ranges, and given that these three ranges sum to ~ 1/3 of 
the total, we take the apparent excess of 30 total counts to correspond to a limit three times 
higher, or ~ 100 counts, on any truly diffuse emission due to the Bondi flow. This limit 
corresponds to an emitted flux of less than 9 x 10~ 15 erg cm -2 s _1 (0.3-7.0 keV) assuming a 
temperature of 0.3 keV, or a luminosity at 780 kpc of < 6 x 10 35 erg s _1 . 

Given that we do not detect a significant enhancement in the flux at (or inside) the 
Bondi radius, we would not expect to detect a temperature change either as any such change 
would need to be co-incidentally offset by a change in electron density. None the less, we 
did select energy bands of Figure 7 such that approximately equal counts are produced in 
the soft and medium bands in order that any change in the ~ 0.3 keV temperature of the 
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Table 2: This table shows the Observation ID's, Dates, Start Times and Exposure Times for 
the 29 ACIS-I images merged together and shown in Figure 7 (left). Total exposure time is 
148 ks. 



diffuse gas would be emphasized. As evidenced by the lack of a color change in Figure 7, 
there is no apparent change in the temperature of the diffuse gas at the Bondi radius. 
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Fig. 7. — (left) The smoothed, summed, 148 ks exposure ACIS-I image centered on M31*. 
The circle is the 5" Bondi radius. The color coding shows the soft band in red (0.3-1.0 keV), 
the medium band in green (1.0-2.5 keV) and the hard band in blue (2.5-7.0 keV). There is 
no clear change in the temperature of the diffuse gas at or inside of the Bondi radius. The 
low temperature of the super-soft source immediately to the South of M31*is apparent by its 
yellowish color, (right): The longest single ACIS-S exposure of M31* is this 40ks exposure 
(obsid 1575). The HST/ACS F435W contours are overlaid and the 5" Bondi radius is shown. 
Pixel size has been set to 0.125" to match the HRC images shown earlier. 

Because spectral fitting summed ACIS data spanning many years consisting of 29 sepa- 
rate observations is not straightforward, we carried out a spectral analysis with the longest 
single ACIS observation of M31*, OBSID 1575. This 40 ks ACIS-S exposure is shown in 
Figure 7 (right), with Bondi radius and also the HST contours overlaid. Extracting the 
counts in a 0.4" radius centered at M31* and excluding the overlapping region centered at 
PI we find 70 counts and find a good fit to a power-law spectrum of the form F V}X oc i>~ a " ,x 
with energy index a UyX = 0.9 ± 0.4 and Nh consistent with Galactic value of 6 x 10 20 cm -2 . 
Freezing Nh to be equal to the Galactic value we are able to place a tighter constraint on 
a VyX = 0.9 ± 0.2. 

Extracting the diffuse counts within the Bondi radius, and excluding the point sources, 
we find 182 counts. Fits to a single component model do not yield acceptable y 2 , but fits 
to a power-law (which is suitable to represent the undetected point sources and the small 
amount of scattering from the detected point sources) and a thermal spectrum are acceptable. 
Freezing the power-law to a U}X = 0.7 and Nh to 6 x 10 20 cm -2 (consistent with the Galactic 
value) yields a good fit (% 2 = 1.3 with 9 dof) with kT= 0.34 ±0.05. The density within this 
region, assuming a 5" radius sphere, is measured to be ~ 0.1 cm~ 3 . 
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2.3. VLA 

We observed M31* simultaneously with the VLA at 5 GHz during our 2004/2005 Chan- 
dra/HRC observations. Figure 8 (right) shows the day-long averages of the radio flux ob- 
served during these observations. We did not see any significant variability within the day 
long averages, to a limit of 15% on timescales of 2, 4, and 6 hours. On other occasions 
M31* has shown variability on several hour long timescales, as can be seen in Figure 8 (left). 
Here the hypothesis of a constant source can be rejected at the 99.99% confidence level 
{x 2 / u = 43.1/11). This high confidence level is due almost entirely to the two measure- 
ments of zero flux in 2002 July 06, as if we disregard those measurements the hypothesis of 
a constant source is acceptable at the 24% confidence level (x 2 ' jv = 11.6/9). 

M31* has been observed at several radio wavelengths from the VLA, but not at the 
same time. Because the source is variable determining a spectrum from these observations is 
uncertain, but it is all the current data will allow. The mean flux densities observed at the 
VLA are ~ 30 ^Jy at 8.4 GHz (measured on 1990 July, 1992 November, 1994 July, 1995 July, 
1996 January), ~ 50 /iJy at 4.9 GHz (measured on 2002 July & August, 2003 June, July 
& August, and during our monitoring simultaneous with Chandra in 2004 December, 2005 
January & February) and between 100 and 140 /iJy at 1.4 GHz (measured on 1981 August, 
1986 August & September). The beam size varies from 0.24" to 0.4" to 1.4", respectively 
at these frequencies, so one must also beware that if there is any diffuse emission it would 
elevate the fluxes at the lower frequencies, however the available images do not show any 
such emission. Taken at face value these flux densities indicate a slope a u ^ r = 0.8, where 
F u , r oc v~ av ' r . 



3. Discussion 

3.1. A Random Superposition? 

While we have detected X-ray emission at a location consistent with M31* which we 
argue does indeed originate from M31*, we must also ask what the chances are of a random 
X-ray binary within the nuclear region of M31 being located at the position of P3=M31*? 
In order to answer this question one needs to estimate the density of sources in the nuclear 
region. As the density increases with decreasing radius one needs to assume some fiducial 
radius within which to compute the density. 

We take the Bondi radius as this fiducial radius. We then count the number of sources 
within this radius both in M31, and also in our Galaxy under the assumption that the density 
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Fig. 8. — (left) The 2-hour averaged radio variability of the nuclear source M31* in three 8- 
hour 6 cm VLA B-configuration observations during 2002 (Sjouwerman et al. 2005). During 
the first half of the second observation M31* was not detected, while it was detected during 
the second half, implying significant variability on a few hour long time scale, (right) The 
day-long VLA averages at 5 GHz from our simultaneous radio/X-ray observations during 
2004/2005. Variability on the few hour long time scale seen in 2002 is not evident in any 
of these 4 individual observations, suggesting that M31*may have been in a different (non- 
variable?) state during 2004/2005. 



of bright nuclear sources is similar. Within the Galaxy, the number of > 10 36 erg s _1 at the 
M31* Bondi radius if it was transferred to the Galaxy (ie, within 100 x 5") is two (Wijnands 
et al 2008). Given this number, the probability of a random source being within 0.2" (the 
registration/centroiding error) of M31* is 0.3%. If we count sources > 10 36 erg s" 1 within 
the Bondi radius at M31 itself, we find 4 (or 5, depending upon the epoch) and a probability 
of 0.6%. We note that these numbers are conservative, because if we increase the fiducial 
radius by 3.2x and therefore the surface area by lOx, the average source density and therefore 
probability decreases by 5x. 

While the odds of an interloper are small, we note that there is nothing in the X- 
ray spectrum of light curve of the source at P3 to distinguish it from an interloping X-ray 
binary. Our identification of this source as M31* is based solely on positional coincidence. 
The variations in luminosity that we have found of 0.6 — 20 x 10 36 erg s " 1 (0. 1-7.0 keV) are 



consistent with the upper limit of 1.2 x 10 erg s set by iLi et al.l (120091 ). particularly 



considering that M31* was relatively faint in the 2004/2005 observations. 



As well as confirming our tentative detection of M31* ( iGarcia et al.l 120051 ) . we have 
detected variable emission consistent with a location in PI. The emission could be due to a 
single low mass X-ray binary embedded in PI. 
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3.2. Bondi Rate 



We have previously estimated the Bondi accretion rate (IGarcia et al.ll2005l ). Herein we 
update that estimation using the same methods. We measure the temperature of the diffuse 
gas within 5" of M31* to b e 0.34 ± 0.05 keV, consistent with what has been previously 
measured in the inner ~ 1' ( ITakahashi et all 12004 ; iLi fc Wans l2007al ; iBogdan fc Gilfanov 



20081 ). This temperature and the new mass estimate lead to a Bond i radius of 5.2" at 
780 kpc. The density in the inner 1" has been estimat ed at 0.06 cm -3 (ILi fc Wangi l2007al ) 
and 0.1 cm -3 ( Takahashi et all2004 ; Dosaj et al. 2001 ). Herein we find a density of 0.1 cm -3 
within the Bondi radius, leading to a Bondi accretion rate of 7 x 10~ 5 M yr _1 and a Bondi 
luminosity of 4 x 10 41 erg s _1 . Given the observed luminosity of 
to an under luminosity of 5 x 10~ 6 as shown in Figure 1. 



2 x 10 3b erg s this leads 



While the Bondi rate is the standard with which to compare SMBHs and to search for 
evidence of RIAFs, it assumes that the gas is stationary at the Bondi radius. Winds from 
stars, su pernovas, or ot her sources may mo d ify th e Bondi rate if they are at high enough 
velocity (IMelial Il992ah . IBogdan fc Gilfanov! (120081 ) suggest that SNR drive a bipolar wind 
out of the plane of M31 with a speed of ~ 60 km s _1 . This wind is seen at large distances 
from M31*, but if it extends to, or originates in, the nuclear region of M31 it still should 
not greatly effect the Bondi flow because its velocity is a factor of five lower than the sound 
speed of ~ 300 km s _1 . On the other hand, if the bipolar flow centers on M31* then there 
may be some connection with outflows or jets originating from the RIAF. 



3.3. Variability 



In the X-ray, the most rapid variability which is clearly detected is a factor of 3 on a 
10 day timescale. This is consistent with the orbital timescale at ~ 100 Schwarzchild radii 
(r$ = 2GM/c 2 ) from M31*. X -ray flares are seen in Sgr A* with durati ons of hours and 
amplitudes of a factor of ~10 (IBaganoff et al.ll200ll ; iMarrone et all 120081 ) . Given that we 
would expect the time scale of such flares to scale with mass, similar flares in M31* would 
occur on a ~ 40 hour or a few day long time scales. Thus this most rapid variability we see 
from M31* might plausibly be associated with the flares seen in Sgr A*. 

Figure 6 shows an (unsuccessful) attempt to detect X-ray variability on even shorter 
time scales. If we had detected variability on very short time scales it might be taken as 
evidence that the X-ray emission was not from M31*, but from an interloping X-ray binary. 
However, we caution here that the 2r$ light crossing time for M31* is ~ 1/3 hour, so one 
would have to detect variations on time scales shorter than this in order to argue against an 
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origin in M31" 



In the radio, we detect more rapid variability than in the X-ray. In particular, during 
the first half of an 8 hour observation on 2002 July 06 M31* is undetected and during the 
second half it is clearly detected. This implies variability of at least a factor of 2 on several 
hour time scales. Sgr A* typically shows onl y 10% to 40% variab il ity in the radio on time 
scales as short as a day iBrown fc Lol (119821 ); iMacquart fc Bowerl (120061 ). Given this, the 
variations we see on 06 July 2002 are unusual. 

The spectral slope of the emission from M31* appears to be the same in the radio and 
X-ray frequency ranges, but the X-ray flux is far above the extension of the radio spectrum 
into the X-ray range. As is the case with Sgr A*, this energy distribution is consistent 
with a synchrotron source for the radio emission and Compton up-scattering of the radio 
photons by t he relativistic electrons i nto the X-ray range in synchrotron self C ompton or 
SSC process jFalcke fc Markoffj I2OO0I : iLiu fc MeliJ l200ll ; Irlornstein et alil2007h . If this is 
indeed the emission mechanism, then the short timescale radio variability should be echoed 
in the X-ray flux. Despite the fact that we did not detect X-ray variability on this short 
timescale during the observations presented herein, the X-ray counting rate when M31* is 
at its brightest is sufficient to suggest that future observations may detect this. 



3.4. Accretion Flow 

Sgr A* appears to be a slightly resolved (non point) source in the Chandra images 
with a luminosity of ~ 10 34 erg s -1 . This implies that the accretion flow itself is bright 
in X-rays. In contrast, in M31* there is no discernible difference in the X-ray emission 
within the Bondi radius. It appears consistent with that from the diffuse gas in the larger 
surrounding area. Unfortunately, it is unclear whether we would even detect the accretion 
flow around Sgr A* if it was at the distance of M31*. The fractional Eddington luminosity of 
a n ADAF flow scales with f ractio nal Eddington mass accretion rate (see for example Figure 7 



of iNarayan &: McClintockl (120081 ). M31* and Sgr A* are accreting at approximately the same 
Eddington scaled rate, so the Sgr A* flow would be be ~ 50x brighter around M31* due to 
its higher mass, or ~ 5 x 10 35 erg s _1 . This is on the order of the upper limit to excess 
diffuse emission due to the Bondi flow that we derived above. However, it is important to 
note that any prediction based on a Sgr A*-like ADAF is uncertain by at least an order of 
magnitude because the 0.1 — 1.0 keV X-ray luminosity of the Sgr A* flow is hidden behind 
10 23 cm -2 of absorption. 



If we accept that the A-stars in P3 are from a recent star formation episode then age 
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of the P3 star cluster may be ~ 200 Myr (IBender et al.ll2005l ). and the mass required to 
produce this cluster during a single star forma tion burst may be between a few 10 4 M and 
10 6 M ( IBender et al.ll2005l ; IChang et al.ll2007l ) respectively. This implies there was a source 
providing this mass at a rate of ~ 10~ 4 to 10~ 2 M yr -1 . We note that the Bondi accretion 
rate is comparable to the lower range of the rate needed to form the star cluster, suggesting 
that the SMBH itself may i nfluen ce the rate of star formation in its immediate surroundings 
(see also iFatuzzo fc Melial (120091 ) concerning this issue with regard to the Galactic center). 
Of course, suggesting that Bondi accretion is important begs the question of where the gas 
that is being accreted came from in the first place. We also note the co-incidence between 
the suggested age of the P3 star clu ster (200 My r ) and the time of the last crossing of M32 
through the nuclear region of M31 iBlock et al.l (120061 ) . suggesting that this crossing may 
have provided the trigger for the star formation event which formed the UV bright cluster 
at P3. 



Conclusions 



We have detected X-ray emission from both optical nuclei of M31 (PI and P3) and 
this variable emission is consistent with two point sources. Having resolved M31* from the 
surrounding point and diffuse X-ray sources we can now investigate the accretion properties 
of this nearby SMBH. The presence of a hot and tr uly diffuse emission compon ent in the 
core of M31 was first note d in Einstein observ ations (jTrinchieri fe Fabbianolll9911) and later 
confirmed with ROSAT JPrimini et al.l ll993h. XMM- Newton Jshirev et al.l l2001ah . and 
Chandra (IDosaj et al.l 120021 ) observations. Any of this gas within the Bondi radius of the 
SMBH will accrete and possibly generate accretion luminosity. In order to compute the 
Bondi accretion rate we use the X-ray observations to estimate the temperature and density 
of this gas. The resulting rate accretion rate would produce a luminosity of 4 x 10 41 erg s _1 if 
the gas radiated with the canonical ~ 10% efficiency Given the observed luminosity of 
~ 2 x 10 36 erg s -1 , M31* is one of the most under-luminous SMBHs known. 

Because M31* has the most highly resolved Bondi flow of any SMBH, a long X-ray 
observation could determine the applicability of ADAF, CDAF, ADIOS, or other models to 
quiescent SMBH accretion. This in turn will tell us the form that black hole accretion takes 
for the vast majority of cosmic time. For example, a 400 ks ACIS-S observation would yield 
2500 counts within the Bondi radius (nearly 10 per pixel), sufficient to divide the region into 
5 annuli and determine accurate temperatures in 5 radial rings. This would be sufficient to 
determine the run of temperature with radius in the RIAF. The various RIAF models all 
predict T (l/r) a , so these data could determine a and therefore the structure of the RIAF. 
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